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ABSTRACT 

The notable absence of radio pulsars having measured magnetic dipole surface held strengths above 
Bo ~ 3 x 10 13 Gauss naturally raises the question of whether this forms an upper limit to pulsar 
magnetization. Recently there has been increasing evidence that neutron stars possessing higher dipole 
spin-down fields do in fact exist, including a growing list of anomalous X-ray pulsars (AXPs) with long 
periods and spinning down with high period derivatives, implying surface fields of 10 14 - 10 15 Gauss. 
Furthermore, the recently reported X-ray period and period derivative for the Soft Gamma-ray Repeater 
(SGR) source SGR1806-20 suggest a surface field around 10 15 Gauss. None of these high-field pulsars 
have yet been detected as radio pulsars. We propose that high- field pulsars should be radio-quiet because 
electron-positron pair production in their magnetospheres, thought to be essential for radio emission, is 
efficiently suppressed in ultra-strong fields (Bo ;> 4 x 10 13 Gauss) by the action of photon splitting, a 
quantum electrodynamical process in which a photon splits into two. Our computed radio quiescence 
boundary in the radio pulsar P — P diagram, where photon splitting overtakes pair creation, is located 
just above the boundary of the known radio pulsar population, neatly dividing them from the AXPs. 
We thus identify a physical mechanism that defines a new class of high-field radio-quiet neutron stars 
that should be detectable by their pulsed emission at X-ray and perhaps 7-ray energies. 

Subject headings: pulsars: general — stars: neutron — magnetic fields — gamma rays: theory — 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

The magnetic fields of pulsars have traditionally been 
inferred through measurement of both a period P and 
a period derivative P . Presuming that the increase in 
period results from spin-down of the neutron star due to 
electromagnetic dipole radiation energy loss, the derived 
surface field is (Shapiro & Teukolsky 1983; Usov & Mel- 
rose 1995) B = 6.4 x 10 19 (PP) 1 / 2 Gauss; this estimate 
assumes a braking index of 3, and a uniformly magnetized 
stellar interior as opposed to the Manchester and Taylor 
(1977) convention of a pure dipole field. A plot of P vs. 
P for radio pulsars with measured P is shown in Figure 0. 
The derived surface fields of these pulsars spans a range 
of five decades from 3 x 10 8 Gauss to 3 x 10 13 Gauss, as 
indicated by the diagonal lines of constant Bq . No radio 
pulsars have inferred fields above this range, even though 
known selection effects for a given pulsar age do not a 
priori prevent their detection (Cordes, private communi- 
cation). Also shown in Figure [j] is the so-called "death 
line," with PP~ 3 = const, (i.e. for fixed open field line 
voltage), to the right of which there are no detected pul- 
sars. The absence of radio pulsars at periods larger than 
5s has been explained as a turn-off in the radio emission 
when the 7-rays emitted by particles accelerated near the 
polar cap cannot produce electron-positron pairs (Sturrock 
1971; Ruderman & Sutherland 1975; Arons 1983). 

In addition to the radio pulsars, we have plotted in Fig- 
ure [l] the AXPs having measured P and P . The AXPs 
are a group of six or seven pulsating X-ray sources with 



periods around 6-12 seconds, which are anomalous in com- 
parison with average characteristics of known accreting 
X-ray pulsars. They are bright, steady X-ray sources hav- 
ing luminosities Lx ~ 10 35 erg s" 1 , they show no sign 
of any companion, are steadily spinning down, and have 
ages t < 10 5 years. Those which have measured P 
(e.g. Mereghetti & Stella 1995; Gotthelf & Vasisht 1998) 
are shown in Table 1 and have derived magnetic fields 
between 10 14 and 10 15 Gauss, assuming that they are 
spinning down due to dipole radiation torques. The soft 
gamma-ray repeaters (SGRs) are another type of high- 
energy source that has recently joined this group of highly 
magnetized neutron stars. SGRs are transient 7-ray burst 
sources that undergo repeated outbursts. There are four 
known SGR sources (one newly discovered: Kouveliotou, 
et al. 1998b), with two (SGRs 1806-20 and 0525-66) as- 
sociated with young ( r < 10 5 yr) supernova remnants. 
Recently, 7.47s and 5.16s pulsations have been discovered 
(Kouveliotou, et al. 1998a, c; Hurley et al. 1998) in the 
quiescent X-ray emission of SGRs 1806-20 and 1900+14, 
respectively, with SGR 1900+14 exhibiting a 5.15s period 
in a 7-ray burst (Cline et al. 1998). In the P — P diagram, 
both AXPs and SGRs live in a separate region above the 
detected radio pulsars. Such a class of highly magnetized 
neutron stars, or "magnetars" , was postulated as a model 
for SGRs (Thompson & Duncan 1993) and more recently 
also for the AXPs (Thompson & Duncan 1996). According 
to the model, magnetars are neutron stars born with mil- 
lisecond periods that generate magnetic fields above 10 14 
Gauss by dynamo action due to convective turbulence. 
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A High B Radio Quiescence Boundary 
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Fig. 1. — The conventional depiction of pulsar phase space, the 
P - P diagram, with filled circles denoting the pulsars listed in the 
Prin r ctr.n Pulcar f!atalr.cmn (Taylnr Manchester & Lyne 1995; see 



also fittp: //pulsar .princeton.edu/), together with their choice of 
a fiducial position for the long period death line. Three gamma-ray 
pulsars arc highlighted, labelled as in the inset. The heavy curve 
denotes the approximate boundary of radio quiescence for photons 
initially propagating along the field, below which pulsars can be ra- 
dio loud, and above which photon splitting suppresses pair creation 
so that pulsars should be radio quiet. Photons were assumed to em- 
anate from near the stellar surface. Also shown on the upper right 
are the positions of four anomalous X-ray pulsars and SGR 1806-20, 
whose measured P and inferred fields Bo are listed in Table 1. The 
curves marked £ ma x > 10 MeV and > 100 MeV are contours (drawn 
only in the radio quiet region) for the photon splitting escape en- 
ergy (see text); to the right of these curves the magnetosphere is 
transparent to 10 MeV and 100 MeV photons, respectively. 

This paper argues that there is a simple physical reason 
(within the confines of polar cap surface emission mod- 
els) why radio emission should be absent in pulsars with 
ultra-strong fields. The behaviour of quantum mechani- 
cal radiation processes in magnetar-type fields is distinctly 
different from that of conventional radio pulsars, where 
Bo <; 10 13 Gauss. While the purely quantum process of 
single-photon pair production 7 — > e + e~ is commonly 
invoked in radio pulsar models (following Sturrock 1971) 
of sources such as the Crab and Vela, at the high field 
strengths appearing in magnetars, the more exotic pro- 
cess of photon splitting, 7 — > 77 , can become quite prob- 
able. The splitting of one photon into two lower energy 
ones (the rate for which was first calculated correctly by 
Adler et al. 1970; Bialynicka-Birula & Bialynicki-Birula 
1970; and Adler 1971), which is third-order (in the fine 
structure constant) in quantum electrodynamics (QED), 
can act as a competitor to pair creation 7 — > e + e~ as 



an attenuation mechanism for gamma-rays. Photon split- 
ting is also a purely quantum effect, and has appreciable 
reaction rates (Adler 1971; Baring & Harding 1997) only 
when the magnetic field is at least a significant fraction 
of the quantum critical field B CI = 4.413 x 10 13 Gauss; 
it forms the centerpiece of the discussion here due to its 
potential for suppressing the creation of pairs. Presuming 
that a plentiful supply of pairs is a prerequisite for, and 
maybe also guarantees, coherent radio emission at observ- 
able levels (a premise of standard polar cap models for 
radio pulsars: e.g. Sturrock 1971; Ruderman & Suther- 
land 1975; Arons 1983), an immediate consequence of the 
suppression of pair creation in pulsars is that detectable 
radio fluxes should be strongly inhibited. We propose here 
that when splitting 7 — * 77 dominates the attenuation of 
gamma-rays, pulsars should be radio quiet. 

TABLE 1 

Spin-Down Parameters' 1 for Anomalous X-ray Pulsars and SGR 1806-20 
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(sec) 


(sec sec -1 ) 
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(Gauss) 
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G10.0-0.3 


7.47 


8.3 x 10"" 


1.4 X 10 3 


1.6 X 10 15 




IE 1841-045^ 
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11.77 


4.7 x lO" 11 


4.0 x 10 3 


1.5 X 10" 




IE 1048-5937 6 




6.45 


2.2 x 10-" 


4.6 x 10 3 


7.6 x 10" 
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8.69 


2.3 x lO" 12 


6.0 x 10" 


2.9 x 10" 




IE 2259+586" 


CTB 109 


6.98 


7.3 x 10- 13 


1.5 x 10 s 


1.4 x 10" 



Note. — (a) Only sources with measured P are included in the Table, (b) The 
magnetic field estimates are obtained using the choice for the dipole spin-down energy 
loss rate B = 6.4 x 10 19 (P P) 1 ' 2 Gauss adopted by Usov & Melrose (1995), derived 
in, foi example, Shapiro & Teukolsky (1983). References for observed periods and 
period derivatives are (c) Kouveliotou et al. (1998), (d) Vasisht and Gotthelf (1997), 
(e) Oosterbroek et al. (1998; here the estimate obtained from their Table 3 using the 
maximum time baseline is used), (f) and (g) Mereghetti & Stella (1995). 



2. CRITERION FOR RADIO QUIESCENCE 

The determination of the pulsar parameter space (i.e. 
Bq , P , P) when such a suppression of pairs by pho- 
ton splitting becomes efficient is therefore of paramount 
interest. To do this, one must compare the attenuation 
probabilities of photons by the two quantum processes, 
7 — > e + e _ and 7 — > 77 , during their passage through the 
pulsar magnetosphere from emission points near the stel- 
lar surface. Polar cap pulsar models invoke an electric field 
parallel to B that rapidly accelerates electrons and/or 
positrons to Lorentz factors of around 10 5 - 10 7 (e.g. Ru- 
derman & Sutherland 1975; Arons 1983; Usov & Melrose 
1995; Harding & Muslimov 1998). If the polar field B is 
not too high (Harding & Muslimov 1998), these energetic 
particles radiate curvature photons (Daugherty & Hard- 
ing 1982), while resonant Compton scattering of quasi- 
thermal X-rays from the stellar surface is the dominant 
radiation production mechanism (Sturner 1995; Harding 
& Muslimov 1998) when B £ 10 13 Gauss. These emis- 
sion processes beam the primary radiation into a very nar- 
row cone virtually aligned with the local magnetic field. 
Non-zero attenuation coefficients T(lo, ^b) (defined as 
the rate divided by c) for 7 — > e + e~ and 7 — > 77 are 
only possible for oblique angles #kB of photon propaga- 
tion with respect to the field: this can be deduced from 
the Lorentz transformation property (Baring & Harding 
1998) T(uj, 6> kB ) =sm0 kB T(wsin0kB, tt/2) for photons of 
energy Lom e c 2 . Hence, to obtain significant optical depths 
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to attenuation, propagation of radiation through the mag- 
netosphere is essential to effect crossing of field lines. 

Both sin 0kB and u> sin #kB generally increase monoton- 
ically during propagation away from the polar cap (Hard- 
ing, Baring & Gonthier 1997). Pair creation cannot take 
place before its threshold of 2m e c 2 / sin#kB is crossed dur- 
ing such propagation. Even though photon splitting has 
a lower rate (being a higher-order QED process), it may 
dominate pair creation because it has no threshold. This 
effect is illustrated in Fig. 1 of Baring & Harding (1997), 
which also exhibits clearly that T(uj sin O^b, tt/2) is a 
strongly increasing function of both u> sin #kB and B for 
both pair creation and photon splitting. 

Computing attenuation probabilities of the two pro- 
cesses for photon propagation near a neutron star requires 
the incorporation of general relativistic effects (Harding, 
Baring & Gonthier 1997; Baring & Harding 1998). Curved 
spacetime modifies photon trajectories from rectilinear 
propagation and therefore changes the values of #kB , red- 
shifts photon energies, and increases the magnetic field in 
the local inertial frame. All three of these effects can lead 
to significant reductions in the attenuation lengths of the 
two quantum processes. Details of such attenuation calcu- 
lations are deferred to Baring & Harding (1998); here we 
focus on their consequences. The first of these is that for 
both processes, due to the r~ 3 falloff of the dipole field, for 
each photon emission point there always exists (Harding, 
Baring & Gonthier 1997) an energy called the escape en- 
ergy, below which the magnetosphere is transparent to the 
radiation, and above which it is opaque and photon escape 
is prohibited. Such domains of opacity at high energies are 
a principal reason why the polar cap model predicts that 
pulsars like the Crab and Vela cannot be visible above a 
maximum energy e max of a few GeV (Daugherty & Hard- 
ing 1982), and perhaps (Harding, Baring & Gonthier 1997) 
also why no emission is observed for PSR1509-58 above a 
few MeV. 

An approximate measure of when the creation of pairs 
is strongly suppressed by photon splitting is when their 
escape energies £ csc are equal. A more accurate measure 
will require detailed simulation of splitting/pair cascades. 
Generally in fields surmised for Crab-like and Vela-like pul- 
sars, e esc for splitting exceeds that for 7 — > e + e~ so that 
photon attenuation only takes place via conversions into 
pairs. In near critical and supercritical fields (B > B CI ), 
the opposite occurs (Harding, Baring & Gonthier 1997; 
Baring & Harding 1998) and splitting assumes the role 
of the dominant attenuation mechanism, strongly (but 
not totally) suppressing the creation of pairs. This ef- 
fect appears to be inevitable, regardless of whether free 
pair creation or positronium formation is considered (Bar- 
ing & Harding 1998), given the polarization properties of 
the photon production mechanisms of curvature and syn- 
chrotron radiation, and resonant Compton scattering, and 
those of the two attenuation processes. 

The condition of equality of escape energies for split- 
ting and pair creation therefore implicitly defines an ap- 
proximate boundary of radio quiescence, which for a given 
colatitude O for the emission of a photon near the sur- 
face, corresponds (Baring & Harding 1998) to a partic- 
ular value of the spin-down dipole field Bq . Such a 
(Bo, O) relationship can be quickly transformed into one 



between the principal radio pulsar observables, namely 
P, the pulsar period (in seconds), and P, the period 
derivative. This is effected using the standard relation- 
ship sinO = [2nR/(Pc)] 1 ^ 2 between the period and the 
polar cap size or colatitude O for a neutron star of ra- 
dius R, modified (Baring & Harding 1998) by geometric 
distortions introduced by curved spacetime, and the dipole 
radiation formula (Shapiro & Teukolsky 1983; Usov & Mel- 
rose 1995) B = 6.4 x 10 19 (PP) 1 / 2 Gauss. The boundary, 
thus transformed, is depicted in Figure [l] in the conven- 
tional P - P diagram for the radio pulsar population, and 
corresponds roughly to 

P « 7.9 x 10-- j (l) 

The particular value —11/15 of the index follows from 
the dependence of the rate of photon splitting on energy 
( tx e 5 ) , photon angles to the field ( oc 9^ B ) , and the mag- 
netic field in this regime (approximately proportional to 

Bq 5 ^ 4 for B cr <^ Bq <^ 4B CI ), coupled with the fact that 
propagation guarantees that #kB oc O , and that the pair 
creation escape energy scales as 9" 1 (see Baring & Hard- 
ing 1998). This putative boundary for radio quiescence, 
which is computed specifically for photon origin near the 
stellar surface, is comfortably located above the entire col- 
lection of radio pulsars, a situation that may be tested by 
the imminent dramatic increase in the population due to 
the new multi-beam Parkes survey (Manchester 1998). A 
study of self-consistent particle acceleration above a pul- 
sar polar cap reveals that emission near the neutron star 
surface is expected for pulsars having B ;> 10 13 Gauss 
(Harding & Muslimov 1998), principally because the pair 
formation fronts are established in regimes where resonant 
Compton scattering dominates the cooling of primary elec- 
trons and created pairs. 

Above the boundary, pair creation is strongly inhibited 
and by extension, pulsars should be radio quiet. The ab- 
sence of radio pulsars in this high- Bq region of phase 
space has become of much greater importance with the 
mounting observational evidence for neutron stars with 
such large fields. This evidence includes the detection of 
spin-down in the growing number of anomalous X-ray pul- 
sars, which when combined with their association with 
young supernova remnants, has shifted the focus from 
accretion torques (e.g. Mereghetti & Stella 1995; van 
Paradijs, Taam & van den Heuvel 1995) to electromag- 
netic dipole torques (e.g. Vasisht & Gotthelf 1997) as the 
origin of the spin-down, implying immense supercritical 
( B > 4.41 x 10 13 Gauss) fields in these sources. These 
are radio quiet pulsars, as is the quiescent counterpart of 
SGR 1806-20, whose 7.47 second periodicity has only been 
very recently announced (Kouveliotou, et al. 1998a). The 
spin-down parameters for these sources, listed in Table 1, 
are used to derive their positions on the P-P diagram 
in Figure |Tj. Our approximate boundary for radio quies- 
cence neatly separates the radio loud and radio quiet pul- 
sar populations. We contend that this pulsar dichotomy 
is a simple consequence of the quantum physics in strong 
fields discussed herein. 
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3. DISCUSSION 

The location of the radio quiescence boundary above the 
detected radio pulsar population gives consistency with 
the assumption of high-energy photon emission at the rim 
of a standard polar cap: larger cap sizes would permit 
photons to acquire large angles to the field more quickly 
pushing the quiescence boundary to lower Bq , in conflict 
with the data. In the extreme case where the emission 
does not arise at the polar caps but occurs in equatorial re- 
gions of the magnetosphere, as advocated (Baring & Hard- 
ing 1995) in SGR models, photon splitting can dominate 
even more efficiently. There are a number of other reasons 
why the emission properties of ultra-magnetized pulsars 
far above the radio quiescence line may well be quite differ- 
ent. The fact that the observed luminosities of AXPs are 
much larger than their spin-down luminosities requires an 
alternative energy source to power their emission. Decay 
of the magnetic field, predicted to be much faster in mag- 
netars than in radio pulsars (Thompson & Duncan 1996), 
may heat the neutron star surface and generate Alfven 
waves that lead to particle acceleration and non-thermal 
emission in the magnetosphere. Since the energy density 
of these fields produces a high level of stress on the neu- 
tron star crust, continual glitches and crustal fracturing 
may occur. The emission from these sources will be gov- 
erned by QED physics of ultra-strong fields, making them 
unique laboratories for the study of exotic processes such 
as photon splitting and even trident production, neutrino 
and muon pair production, or magnetic Cerenkov radia- 
tion, in addition to more commonplace mechanisms such 
as cyclotron/synchrotron radiation and resonant Compton 
scattering. 

The observational expectations for high-field pulsar 
emission in the radio, X-ray and gamma-ray bands can be 
quickly summarized in view of the physical effect of sup- 
pression of pair creation. Clearly, radio luminosity should 
decline as the boundary of radio quiescence is approached 
from below, an effect that can be enhanced by the creation 
of pairs in low Landau levels (Harding & Daugherty 1983) 
or positronium formation (Usov & Melrose 1995) above 
Bq j> 6 x 10 12 Gauss. Note that photo-ionization and 
electric field dissociation of positronium can mute the ef- 
fectiveness of bound state pair creation in inhibiting radio 



emission (such issues are discussed in Baring & Harding 
1998). Since synchrotron radiation, the dominant emis- 
sion of produced pairs, has a steeper spectrum than pri- 
mary curvature or resonant Compton-scattered radiation, 
the absence of a synchrotron component should lead to a 
progressive flattening of pulsar gamma-ray spectra as Bq 
is increased to values above ~ 4 x 10 13 Gauss. In the mag- 
netar regime of ultra-strong fields, the maximum energy 
of emission e max is dictated by the escape energy e csc for 
splitting; sample values are illustrated in Figure [D, indicat- 
ing that SGR1806-20 in its quiescent state and the Kes 73 
and IE 1048-5937 AXPs would barely be visible at around 
100 MeV. A decline in gamma-ray luminosity below the 
Bq I P 2 proportionality found for conventional gamma-ray 
pulsars (e.g. Thompson, ct al. 1997) would also be ex- 
pected. In the X-ray range, non-thermal spectra would 
probably also flatten as Bq increased above 10 13 Gauss, 
as the cyclotron energy moved into the MeV range. How- 
ever, at much higher fields, the AXPs exhibit steep X-ray 
continuum emission (and probably distinct from any hard 
X-ray or gamma-ray component), at luminosities far ex- 
ceeding that available due to dipole spin-down. Hence, 
given that another mechanism for energy supply such as 
field decay (Thompson & Duncan 1996) must be operat- 
ing in these sources, how their spectrum and luminosity in 
the X-ray band would depend on Bq and P is not well 
understood. 

The prediction of a radio quiescence boundary, coupled 
with the expectations mentioned just above for emission 
properties in the soft gamma-ray band, provide strong ob- 
servational tests for the quantum effects in strong mag- 
netic fields relating to pair production and its suppression 
by photon splitting. Future studies by radio telescopes 
and measurements by X-ray experiments such as AXAF 
and XMM, and gamma-ray missions such as Integral and 
GLAST should delve deeper into the exciting phenomenon 
of pulsars with ultra-strong magnetic fields. 

We thank Alex Muslimov, Jim Cordes and Eric Got- 
thclf for discussions. MGB acknowledges the support of a 
Compton Fellowship during the period when part of this 
research was completed. This work was supported by the 
NASA Astrophysics Theory Program. 
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